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D ' We present the thermodynamic properties of strongly interacting matter in finite volume in the 

framework of Polyakov loop enhanced Nambu— Jona-lasinio model within mean field approxima- 
tion. We considered both the 2 fiavor and 2+1 fiavor matter. Our primary observation was a 
\l ' qualitative change in the phase transition properties that resulted in the lowering of the tempera- 

, ture corresponding to the critical end point. This would make it favorable for detection in heavy-ion 

■ experiments that intend to create high density matter with considerably small temperatures. We 

(— I ' further demonstrate the possibility of obtaining chiral symmetry restoration even within the confined 

phase in finite volumes. 

Oh' PACS numbers: 12.38.Aw, 12.38.Mh, 12.39.-x 

CD ■ 

^ ; I. INTRODUCTION 

cn ■ 

The strongly interacting matter is supposed to have a rich phase structure at finite temperatures and densities [l| . 
While our Universe at the present epoch contains a significant fraction of color singlet hadrons, color non-singlet 
states especially quarks and giuons may have been prevalent in the past — a few microseconds after the Big Bang [2] . 
. The temperature of the Universe at that epoch is estimated to be ^ 200 MeV. Similar state of matter is also expected 
to exist inside the core of super-massive stars in the present day Universe, where the density is ~ 10 times that of 
normal nuclear matter. A direct study of such natural phenomenon is out of bounds even to modern astrophysicists. 
Fortunately experimental facilities at CERN (France/Switzerland), BNL (USA) and recently at GSI (Germany) are 
exploring the possibilities of creating and studying the properties of such exotic states of matter in a controlled 
environment. The key differences that appear in such experiments as compared to the natural phenomenon arc the 
lifetime of matter created in the exotic state and its volume. Whereas in natural phenomenon the lifetime of the exotic 
• matter may be large compared to the interaction time-scale, it is usually very small in an experimental situation. 
Some effects of the enhanced lifetime on the physical aspects of the system relating to the onset of equilibrium of weak 
interactions were discussed by us in Ref. 0] . Here we shall discuss about the effects of finite volume on the properties 
of strongly interacting matter. 

In the following we shall generically define the matter with color confined states as the hadronic phase and the exotic 
state with colored degrees of freedom as the quark-gluon plasma (QGP). In the experiments this exotic phase may be 
produced by ultra relativistic collisions of heavy ions. The volume of the system thus created would depend on the 
nature of the colliding nuclei, the center of mass energy (y^) and the centrality of collision. Once created, the system 
expands until the constituents are so far separated that their interaction ceases and they flow out as free streaming 
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particles. The distribution of particles thus freezes out, except for some further decays to smaller particles. There 
have been a large number of efforts to estimate the system size at freeze-out for different yfs and different ccntralities. 
The most popular way of doing so is to measure the HBT radii. In Ref. [3l it has been shown that the freeze out 
volume increases as the -^/i increases. Here the authors have estimated the freeze out volume and found it to vary 
from 2000 fm^ to 3000 fm^. In a very recent paper Q the volume of homogeneity has been calculated using UrQMD 
model ^ and compared with the experimentally available results. The -y/s considered was in the range of 62.4 GeV to 
2760 GeV for lead-lead collisions at different ccntralities. The system volume has been found to vary from 50 frm? to 
250 fnv^ . Given that these are the freeze-out volumes, one can trace back to the initial equilibration time and expect 
an even smaller system size. In fact one cannot even consider the whole fireball, which is an isolated system to be in 
thermodynamic equilibrium. One has to choose a proper rapidity interval to act as the system under consideration. 
Therefore it becomes important to study how the various thermodynamic quantities in a strongly interacting matter 
depend on the volume of the system. Specifically we know that finite system sizes would lead to smoothening of any 
singularities appearing at a phase transition Q. Thus important signatures of such transitions must be reanalyzed 
with the help of finite size scaling analysis Q . In the context of heavy ion collisions such a possible analysis has been 
discussed in the literature (sec e.g. @t11|)- 

On the theoretical side a study of finite volume effects was done in Ref. [l^l with a bag of non- interacting quarks 
and gluons and it was found that the effective degrees of freedom are reduced. In Ref. [13| a two model equation of 
state was used to show that the separation between the hadronic and QGP phases around the critical temperature 
looses its sharpness resulting in a soft effective equation of state. A few first principle study of pure gluon theory on 
space-time lattices were performed, showing the possibility of significant finite size effects [13, [l^. Similar studies 
are going on in various QCD inspired models. In Ref. [la, [l^ the quark mass gap equation has been studied with 
Schwingcr-Dyson equation parallel to equivalent Lattice QCD (LQCD) calculations and various meson properties 
arc found to have significant volume dependence. There are studies with four-fermi type interactions in the NJL 
models [1, [H, [H, linear sigma models [H [13 [U and Gross-Neveu models [22|. While in Ref. ^ the scaling 



behavior of chiral phase transition for finite and infinite volumes has been studied, the character of phase diagram has 
been studied in Ref. [^, [13, [HI, [13 ■ I'^ T:efs. [l^l and [l3| the authors have studied the chiral properties as a function 
of the radius of a finite droplet of quark matter. The stability of such a droplet in the context of strangelet formation 
within the NJL model has been addressed in Ref. [13| • Size dependent effects of difermion states within 2-dimensional 



NJL model has been studied in Ref. [2J| and that of magnetic field is discussed in Ref. [25 1. 

In this work we shall use the Polyakov loop enhanced NJL (PNJL) model to study the thermodynamic properties 
of the strongly interacting matter in a finite volume. This model originated from the NJL model [26l429l | which 
incorporates the global symmetries of QCD quite nicely. A four quark interaction term in the NJL Lagrangian is 
able to generate the physics of spontaneous breaking of chiral symmetry — a property of QCD which is manifested as 
the non-degenerate chiral partners of the low-mass hadrons. However a reasonable description of the physics of color 
confinement is missin g. W ith the introduction of a background field in the NJL model, motivated by the dynamics 
of the Polyakov Loop [30| , one obtains the PNJL model which describes a number of features of confinement physics 
quite satisfactorily (see e.g. (3TI - [36j ). 

Certain aspects of finite volume effects in the PNJL model has been discussed in Ref. [13] through a coarse graining 
of the Lagrangian, followed by a Monte Carlo simulation. This method goes on similar lines as the numerical studies 
of LQCD. Normally this would involve the same kind of complex determinant problem that has plagued the direct 
LQCD computations for non-zero baryon number densities. So it may be desirable to keep using the saddle point 
approximation in PNJL model to study the finite volume effects. Here we make the first case study, albeit with some 
simplified assumptions towards that direction. 

We organize our paper as follows. In the next section we briefly describe the PNJL model and the modifications 
for finite volume. In section IIIII wc describe phase transition at finite volume and in section IIVI we discuss the 
thermodynamic properties. The pion and sigma meson masses and the pion decay constant at finite volume have 
been discussed in section fVl In section IVTl we summarize and conclude. 

II. THE PNJL MODEL 

We shall consider the PNJL model with light fiavors (2 fiavor) and light plus strange flavors (2+1 flavor). In the 
PNJL model the gluon physics comes into play through the chiral point couplings between quarks (present in the NJL 
part) and a background field which represents Polyakov Loop dynamics. The Polyakov line is represented as, 

L{x) = Pexp[i dTA4{x,T)] (1) 
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where A4 = iAo is the temporal component of Eucledian gauge field {A, A4), /? = y , and V denotes path ordering. L{x) 
transforms as a field with charge one under global Z(3) symmetry. The Polyakov loop is then given by $ = {TrcL)/Nc, 
and its conjugate by, $ = {TrcL'^)/Nc. The gluon dynamics can be described as an effective theory of the Polyakov 
loops. The Polyakov loop potential can be expressed as, 

u'imMAiT) ^ uimM (2) 

where L{{(f>) is a Landau-Ginsburg type potential commensurate with the Z(3) global symmetry. Here we choose a 
form given in Ref. [3^ . 

where 

b,{T) = ao + ai(^) + a,i^r + asi^f, (4) 

&3 and &4 being constants. The second term in Eqn.Q is the Vandermonde term which replicates the effect of SU(3) 
Haar measure and is given by, 

J[$, $] = (27/247r2) [f - 6$$ + 4($^ + l>^) - 3($$)^ 

The corresponding parameters were earlier obtained in the above mentioned literature by choosing suitable values by 
fitting a few physical quantities as function of temperature obtained in LQCD computations. The set of values chosen 
here are, 

ao = 6.75, ai = -1.95, 02 = 2.625, 03 = -7.44, 63 = 0.75, 64 = 7.5, 
To = 190 MeV, k = 0.2 (for 2 flavor), k = 0.13 (for 2+1 flavor) 

For the quarks we shall use the usual form of the NJL model except for the substitution of a covariant derivative 
containing a background temporal gauge held. Thus the 2 flavor version of PN JL model is described by the Lagrangian, 

f=u,d f / 

+ f E [{i'r'^^f + {^n5r''^f]-^'{m],m],T) (5) 

a=l,2,3 

For 2+1 flavor the Lagrangian may be written as, 

^= E '/^/7p*^''^/-E"V^/V^/+Ew7oVi/V'/ + Y E [(^A»^ + (^Z75A''^)2] 

f=u,d,s f f a=0,...,8 

- goidet^fPLiJf' + deti^fPnijf,] - W'($[^], T) (6) 

where / denotes the flavors u or d or s respectively. The matrices Pt.fl — (1 i 75)/2 are respectively the left- 
handed and right-handed chiral projectors, and the other terms have their usual meaning, described in details in 
refs. [H, [13, S, I38l - l40l |. This NJL part of the theory is analogous to the BCS theory of superconductor, where the 
pairing of two electrons leads to the condensation causing a gap in the energy spectrum. Similarly in the chiral limit, 
NJL model exhibits dynamical breaking of SU{Nf)j^ x SU{Nf)j^ symmetry to SU{Nf)v symmetry {Nf being the 
number of flavors). As a result the composite operators ipf^pf pick up nonzero vacuum expectation values. The quark 
condensate is given as, 

(ipfijf) -iNcCty^^+{trSf{x - y)), (7) 
where trace is over color and spin states. The self-consistent gap equation for the constituent quark masses are. 



ruf - 55(7/ + 5L10-/+10-/+2, 



(8) 
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where ct/ = {'ipfipf) denotes chiral condensate of the quark with flavor /. Here if we consider cr/ = cr„, then (T/+i = aj, 
and IT/+2 = fs- Similarly if af = ad then 0"/+! = as and cr/_|_2 = (^u, if o"/ = ag then trj+i = cr„ and cr/+2 = ft;. The 
expression for ct/ at zero temperature (T = 0) and chemical potential {fif — 0) may be written as j35| . 

A being the three-momentum cut-off. This cut-off have been used to regulate the model because it contains dimen- 
sionful couplings rendering the model to be non-rcnormalizable. 

Due to the dynamical breaking of chiral symmetry, N'j — 1 Goldstone bosons appear. These are the pions and kaons 
whose masses, decay widths etc. from experimental observations are utilized to fix the NJL model parameters. The 
parameter values have been listed in table HI Here we consider the $, $ and af fields in the mean field approximation 
(MFA) where the mean field arc obtained by simultaneously solving the respective saddle point equations. 



Model 


niu 




A 


gsA^ 


9dA^ 




MeV 


MeV 


MeV 






2 flavor 


5.5 





651 


4.27 





2+1 flavor 


5.5 


134.76 


631 


3.67 


9.33 



TABLE I: Parameters of the Fermionic part of the model. 



Now that the PNJL model is described for infinite volumes we discuss how we implement the finite volume con- 
straints. Ideally one should choose the proper boundary conditions — periodic for bosons and anti-periodic for 
fermions. This would lead to a infinite sum over discrete momentum values pi — nni/R, where i = x,y, z and rii are 
all positive integers and R is the lateral size of a cubic volume. This implies a lower momentum cut-off Pmin = tt/ R = A 
(say). One should also incorporate proper effects of surface and curvatures. In this first case study we shall however 
take up a number of simplifications listed below: 

• We shall neglect surface and curvature effects. 

• The infinite sum will be considered as an integration over a continuous variation of momentum albeit with the 
lower cut-off. 

• We shall not use any modifications to the mean-field parameters due to finite size effects. Thus the Polyakov 
loop potential as well as the mean-field part of the NJL model would remain unchanged. They shall feel the 
effect of changing volume only implicitly through the saddle point equations. 



III. PHASE TRANSITION 



To study the finite volume effects on the thermodynamic properties of strongly interacting matter we begin by 
writing down the thermodynamic potential in MFA. 
The expression is given by, 

n'{^,^,af,T,^lf)=U'[^,^,T] + 2gs ^} ^ ^'^u'Jd'Js 

where a;„ = TrT{2n + 1) are Matsubara frequencies for fermions. The inverse quark propagator is given in momentum 
space by 

S-^ ^joip° + fi-iAi)-^-p-M (11) 
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FIG. 1: (Color online) Derivatives of order parameters for chiral and deconfinement phase transition for different system sizes. 



using the identity Tr In {X) = Indet {X), we get, 

^ . ^ 4^ J\ 27r 



f—u^d^s 



/ 



In 



1 + 3(ci> + $ exp(^^%l^)) exp(^^%^) + cM~^^\~ 



T 



T 



1 + 3($ + $exp( 



T 



-))exp(- 



T 



) + exp(- 



T 



= n- kT^ In J[$,$] 



(12) 



where E, 



Pf 



Mj is the single quasiparticle energy. In the last line contains all the terms of O' except the 



Vandermonde term. 

We now search for the saddle point of the thermodynamic potential which gives the temperature and density 
dependence of the fields. For all the system sizes, at zero baryon density, we found that the order parameters for 
both chiral (cr uu > + < dd >) and deconfinement ($) transition smoothly passes from the hadronic phase to the 
quark phase. This indicates that the system does not have a real phase transition, rather there is a smooth crossover. 
The crossover temperature is identified to the point of inflection of and $ with temperature. In Fig. [T] we have 
plotted d^/dT and dcju/dT for 2 flavor and 2+1 flavor matter for different system sizes. The peak position of these 
plots give respective inflection points. Note that, the deconfinement and chiral transitions do not take place exactly 
at the same temperature. Here we take the average of these two temperatures as Tc. The results are shown in table 
im where we quote the different values of the crossover temperatures corresponding to different system sizes. 

From table HI] it can be seen that the Tc has a strong dependence on system size. For 2 flavors the Tc varies from 
167 MeV to 186 MeV which means a change of about 10%. A similar result is observed for 2+1 flavor. One should 
note that the shift in the Tc is mainly due to the shift in the transition temperature of the chiral phase transition. 
The transition temperature of the deconfining phase transition almost does not change. This result is similar to that 
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R 


= 2 /m 


R 


= 2.5 fm 


R 


= ?> fm 


R 


^5 fm 


R = oo 


Tc (MeV) (2 flavor) 




167 




171 




180 




184 


186 


Tc (MeV) (2+1 flavor) 




160 




167 




174 




180 


181 



TABLE II: Transition temperatures for different system sizes. 



obtained with PNJL model on the lattice [37|. This is somewhat expected as the Polyakov loop potential feels the 
effect of changing volume only indirectly through the fields $ and $. 





FIG. 2: (Color online) Constituent masses of quarks as a function of temperature for different system sizes. 

In Fig. [2] we have plotted the temperature dependence of the constituent quark masses for both 2 flavors and 2+1 
flavors. Below the crossover temperature they exhibit very strong volume dependence. Smaller the volume, smaller 
is the constituent mass. In the 2+1 flavor case, the masses of the light flavors drop faster than the strange quark. 
It thus seems that the chiral symmetry is gradually getting restored as one looks into smaller and smaller volumes. 
This is also the reason why the Tc itself is lowered for smaller volumes given in table |lll Similar feature has also been 
observed in NJL models [isl . [l9j . Given that the quark condensation is similar to the superconducting condensate it 
is interestin g to note that there are in fact certain superconductors which show similar decrease of band gap with the 
system size pTj . 

Let us now take a look into the situation at non-zero quark chemical potential fiq = fXf /Nj. For infinite volume 
the phase transition is of first order and one observes a gap in the order parameter at sufficiently high chemical 
potential. At some smaller /ig, the first order transition ends at a critical end point (CEP). At this point the system 
undergoes a second order transition. At even smaller fiq we have only a crossover. As the volume of the system is 
lowered we find the phase transition characteristics fade away. Even the crossover characteristics start to die down. 
This is clear from the Fig. [3] where we plot duu/dT and d^/dT for /.(^ = 300 MeV as a function of temperature. In 
Fig. |4]the phase diagram as a function of system size is shown. Note that the CEP gradually shifts towards higher 
fig and lower T and finally disappears as the volume is reduced. This is an encouraging fact for the critical point 
search in heavy-ion collision experiments. To attain such high densities one needs to collide the ions at low y/s, which 
means the temperature attained is lower. So if it were an infinite system one would have been far away from the CEP. 



7 




FIG. 3: (Color online) Derivatives of order parameters for chiral and deconfinement phase transition for 2 flavor at finite fj.q. 



Fortunately the experiments would produce small system volumes and this may lead to the location of the respective 
CEP possible. Thereafter one would need to do the finite size scaling analysis to extrapolate to the CEP for infinite 
volumes. The location of CEP for different volumes is collected in table lllll 
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FIG. 4: (Color online) Phase diagram for different system sizes. The inside curve is for R = 2 fm, the next curve is for 
R = 2.5 fm and the outermost curve is for R = oo. 





R = 2 fm 


R = 2.5 fm 


R = 3 fm 


R = 5 fm 


R — oo 


Tc (MeV), fig^ (MeV) (2 flavor) 
Tc (MeV), Atq^ (MeV) (2+1 flavor) 


No CEP 
No CEP 


32, 339 
32, 328 


52, 335 
60, 324 


69, 330 
86, 316 


72, 329 
93, 313 



TABLE III: Location of chiral CEP for different system sizes. 



IV. THERMODYNAMICS 



In this section we discuss the behavior of a few thermodynamic observables namely pressure, energy density, specific 
heat, speed of sound etc. for different system sizes. 

The pressure inside a volume V may be written as, P{T, fiq) = —^i^KZ^alKl ^ where T is the temperature and fXq is 
the quark chemical potential. In the top left panel of Fig. [5] we plot the temperature dependence of scaled pressure 
(P/T-*) for 2 flavor system. As can be seen there is a significant change in scaled pressure for small system sizes. For 
example at Tc the P/T^ for a system with R = 2 fm is almost half of that of an infinite system. As the temperature 
increases the difference slowly diminishes. The decrease of scaled pressure with increasing volume may be a surprise 
given that the constituent quark masses were shown to decrease drastically with decreasing volume, which should 



8 




imply increase in pressure. This can be understood as follows. With decreasing volume, not only the constituent 
masses decrease, but also the lowest momentum increases due to the infrared cut-off. These two conditions somehow 
seem to keep the lowest available energy of the quark quasi-particles almost same for different volumes. Thus the 
pressure does not increase with decreasing volume. However when plotted against T/Tc it seems to decrease because 
the Tc itself is smaller for smaller volumes, and therefore the pressure at the corresponding T/Tc for smaller volume 
is smaller than that for a larger volume. 

ri. In the top 



The volume dependence is also quite strong for the energy density e = — T^^^^^^ 



dT 



right panel of the Fig. [5] we have plotted the e/T^ as a function of T/Tc for different system sizes. It has similar 
characteristics as P/T* but the difference seems to diminish faster with increasing temperature. As the system size 
becomes i? = 5 /m both the scaled pressure and scaled energy density converge to the i? — > oo case for almost all 
temperatures. 

It is well known that for infinite volumes the definition of pressure simplifies to, P{T,^q) = — fl{T,^q), which is 
commonly used the literature for PNJL models at infinite volumes. However since we are considering finite volumes 
here it would be interesting to check how much difference will it make if we keep using this definition rather than the 
correct one with a volume derivative. In the bottom two panels of Fig. [5] we have made a comparison of —fl/T"^ and 
P/T'^. For R = 2 fm we see that these two quantities differ by about 10%. Again, as the size goes close to i? = 5 fm 
this difference is almost washed out. 

Let us now consider the quantity e — 3P. In our mean field approach this is the trace of the energy-momentum 
tensor given by, 7^^ = e — 3P. In a conformaly symmetric theory, for example a theory of free massless quarks and 
gluons the energy momentum tensor is supposed to be zero as it signifies the conservation of the conformal currents. 
Thus e = 3P in that limit. In QCD however the conformal symmetry is broken due to non-zero quark masses as 
well as quantum anomalies as evident from the presence of a scale in the running coupling constant (42l . l43j . Thus 
the energy-momentum tensor does not remain traceless. This was also found to be true in the PNJL model that 
have been reported in our earlier studies and compared with LQCD results [13, The PNJL model is however 
not QCD and the reason for the scale symmetry breaking is the introduction of an ultraviolet cut-off in the NJL 
part, a temperature scale Tq in the Polyakov loop part and of course a quark mass term similar to that in QCD. The 
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FIG. 6: (Color online) Conformal measure C = 1 — 3P/e for different system sizes as a function of temperature. 



physical implication of the two different scales in the quark and Polyakov sector is to give rise to separate crossover 
temperatures for the two sectors. To compare quantities obtained in PNJL model against LQCD results one then 
averages out two crossover temperatures as done by us here in the last section. Now for finite system sizes we have 
introduced an infrared cutoff which should further enhance the effect of conformal symmetry breaking. In fig. [5] we 
show the variation of the conformal measure C = [e — 3P)/e with temperature for both 2 flavor and and 2+1 flavor 
matter for different system sizes. That the smaller system sizes lead to larger conformal symmetry breaking effects 
is evident, except for the anomalous behavior of the lowest size of i? = 2 fm. (Though not shown in the figure we 
found that the anomalous behavior starts at a size between 2 fm < i? < 2.5 fm. The reason for this behavior is not 
clear at the moment and requires further investigation.) 
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FIG. 7: (Color online) Variation of specific heat with temperature for different system sizes. 
The specific heat at constant volume Cy = & 



is shown in Fig. [T] We find that with the change in volume, Cy 

changes prominently up to the temperature corresponding to the crossover region. For smaller volumes the specific 
heat is smaller indicating a higher rise in temperature for the same rise in energy density. Obviously this can be 
correlated with the temperature dependence of energy density discussed earlier in Fig. [51 We found that a given 
amount of scaled energy density is obtained at a higher scaled temperature for a smaller volume. This can be of 
interest in heavy-ion collision experiments. A given energy density deposited in the finite volume would create a 
plasma with temperature somewhat higher than that expected in a similar volume inside an infinite volume system 
having the same energy density. 



The specific heat is also a measure of energy fluctuations in the system [4J] . Fluctuations tend to rise sharply near 
a phase transition. For a crossover they are somewhat subdued. Obviously for finite volumes a true phase transition 
is not possible and as one keeps on decreasing the volume all signatures even for a crossover should die down. This 
is exactly the behavior of Cy as presented in Fig. [71 

The squared speed of sound v"^ = ^ is shown in Fig. [HI At large temperatures the v'^ reaches its maximum value 
as the system becomes almost ideal. Interactions grow with decreasing temperatures resulting in the lowering of v^. 
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FIG. 8: (Color online) Variation of squared speed of sound with temperature for different system sizes. 



The conformal measure C may be considered as a measure of the strength of the interaction in the system. Thus 
lower the value of C, higher should be the value of v1. This is evident from Fig. [HI where we find the v1 to decrease 
with decreasing temperature, just opposite to the behavior of C shown in Fig. [HI This correlation between C and v1 
also apparent for variation in volume. With decreasing volume the speed of sound decreases. (In fact an anomalous 
behavior for the smallest size R = 2 fm is also apparent for uj.) A smaller speed of sound for smaller volumes would 
mean a slower flow for finite size systems created in heavy- ion collisions. 



V. PROPERTIES OF NON-STRANGE MESONS 

For infinite volumes the meson properties in the PNJL model has been discussed for 2 flavors [4^1 as well as for 
2+1 flavors [s^, H^. In this section we describe the properties of non-strange mesons at flnite volumes in the PNJL 
model. A detailed account of the calculational procedure for meson masses at finite temperatures and densities in the 
PNJL model may be found in Ref. (sH]. Here we sketch the outline of the task. 

The collective excitations, the fluctuation of the mean field around the vacuum can be handled within the Random 
Phase Approximation (RPA) [U. In this approximation, which is equivalent to summing over the ring diagrams, the 
retarded correlation function for a meson M is given by, 

^ 1 - 2GmU^' ■ ^^^^ 

Here G]\j is the suitable coupling constant and Hj\/(fc^) is the one- loop polarization function for the mesonic channel 
under consideration. Within the RPA, H^^ may be written as, 

H^^= / ^TT[rMSip + q)TMS{q)], (14) 

where S{p) is the Hartree quark propagator, Tm is the appropriate combination of gamma matrices of different 
mesonic channels and the trace is taken over the Dirac and color spaces. The lower limit on the integration is now 
required for finite volume studies. 

Here we concentrate on the scalar (cr) and pseudoscalar (tt) channels. These contributions can be written as, 



n.(g')= ^Tr {S{p + q)S{q)). (15) 

The pole mass can be obtained by solving, 

l-2GMll"{uj = mM,q^0) = 0. (16) 
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FIG. 9: (Color online) Variation of meson masses (upper panel) and pion decay constant (lower panel) with temperature for 
different system sizes. In the upper panel the three monotonously rising curves are for and the other three are for mo-. 



where niM is the mass of a particular meson. The detailed expression for 11*'' and Gm for tt and a mesons may be 
found in Ref. [36| . 

In the upper panels of Fig. [H]we have plotted the masses of pion (to^) and sigma (too.) as a function of temperature 
for different system sizes. In any given volume we see that for low temperatures the masses of pion and sigma are 
different and they become degenerate above Tc where chiral symmetry is expected to get restored. With decrease in 
volume we find the pion mass to increase. However above 1.2 Tc the pion mass for infinite volume suddenly shoots up 
above those for the finite volumes. This may have important consequences in heavy-ion reactions where system size 
is small. Whereas for infinite volume the fast increasing mass of pion would drastically reduce chances of obtaining 
pion-like bound states, the same may not be true for finite volume. 

While the mass of pion increases with decreasing volume at low temperatures the mass of sigma is found to decrease 
and that too quite fast. One can actually see a trend to the masses of the two chiral partners becoming closer to each 
other with decreasing volume. This, yet again, shows that chiral symmetry breaking effects reduce with decreasing 
volumes. 

The pion decay constant may be obtained from the matrix clement (0| 5|7r''(fc)) — iSabf-jik^^ where J^s = 

ijj^j5-^ip is the chiral current. At finite temperature and for a particular volume it can be written as (see (27l. [28j 
and including the low momentum cut-off A), 



r2 _ 

27r2 



-2 



El Jx El 



(17) 



where Ep = \Jp^ is the single particle energy of a light quark and f{Ep) is the distribution function properly 

modified due to Polyakov loop interaction. 

As shown in the lower panels of Fig. |9l the pion decay constant decreases both with the decrease in temperature and 
with that of system size. The decrease of with temperature has also been observed in other effective models (45l.[48|. 
Dyson-Schwinger approaches [i^ as well as in LQCD [H^l- This is also an indication of the restoration of chiral 
symmetry as is directly proportional to the divergence of the chiral current. 
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The tendency of chiral symmetry getting restored in finite volumes may also be noted by comparing Fig. |9] with 
Fig. O At low temperatures the constituent quark masses decrease with decreasing volume. It so happens that the 
light constituent quark masses become smaller than the pion mass for the smallest sizes studied here. These quarks 
should then become thermodynamically more favored than the pions. Though fortunately in the PNJL model, such 
constituent quarks will be suppressed due to the presence of the Polyakov loop, the pions would still loose their 
significance as the lightest particles that made them suitable candidates for becoming the Goldstone bosons. Thus 
what seems to happen is that the decrease of volume restores the spontaneous breaking of chiral symmetry in the 
same way as increase in temperature. The critical size Rc for such symmetry restoration would be somewhere between 
2 fm and 2.5 fm. From Fig. [2] it may be noted that this range of sizes is almost equal to the respective constituent 
quark masses. This observation is commensurate with the expectation from chiral perturbation theory that chiral 
symmetry restoration may take place once the quark masses become equal to the inverse of system size [5l| . 
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FIG. 10: (Color online) Checking for violation of the GMOR relation as a function of temperature and system size. 



With all the strong indication of a possible chiral symmetry restoration with decreasing volume it would be inter- 
esting to see what happens to the Gell-Mann Oakes Renner (GMOR) [s^l relation, which in the lowest order of chiral 
expansion is given by /^mj = m„(T„ + 771^17^. Normally with increase in temperature as the spontaneously broken part 
of the chiral symmetry gets restored the GMOR relation should start to break down. This is exactly what we find in 
our calculations and shown in Fig. 1101 But surprisingly we find that similar effect is not observed for the decrease in 
volume and the GMOR relation holds good for the all the ranges of volumes we considered almost up to temperatures 
as high as 0.8 Tc. In fact even at higher temperatures the GMOR relation is violated the most for infinite volumes. 
The way one can understand this is that for a physical chiral expansion a quantity m-^/Ai is required, where M is 
some suitable scale. In chiral perturbation theory M is usually the neucleon mass, in zero temperature NJL or PNJL 
models it is the high momentum cut-off A, etc. For finite temperatures one can then consider T to play the role 
of Ai. Thus, given a temperature if the corresponding m^r in a given volume is less than T, chiral identities would 
work properly (see e.g. [53|). So here we have a situation where chiral symmetry is getting restored while partial 
conservation of axial current is still maintained. 



VI. CONCLUSION 



We have tried to understand the dynamics of strongly interacting matter inside finite volume in the framework 
of PNJL model with saddle point approximation. Several interesting results were observed that can have important 
implications for heavy-ion collision experiments. Our major finding was that the spontaneously broken chiral symmetry 
may be restored at much lower temperatures in small volume. This was shown through the computation of various 
thermodynamic observables as well as certain hadron properties. 

Changes in the equation of state and speed of sound may have important consequences in the flow properties of the 
exotic medium created in the experiments. A measure of the specific heat in heavy-ion experiments is the transverse 
momentum fluctuations. We find the specific heat to decrease with decreasing volume indicating that the momentum 
fluctuations may not be as large as expected at a given T/Tc- 

From the variation of the phase boundary with changing volume we demonstrated a stronger possibility of finding 
the signatures of a critical end point in low energy experiments that intend to create high baryonic densities where 
the expected temperature is not too high. 



13 



Finally from the hadron properties we observed the possibility of obtaining a chiral symmetric but confined phase 
in small volumes. We hope that a combination of heavy ion collisions and not-so-heavy ion collisions at similar center 
of mass energies, followed by an appropriate finite size scaling study may give us a better understanding of the QCD 
phase structure. 

As discussed earlier we made a couple of simplified assumptions in this work. The Polyakov loop potential used 
here does not have an explicit volume dependence. The discrete momentum states in the quark potential was replaced 
with a continuum, and the only explicit dependence on system size was through the lower momentum cut-off. Though 
we believe that these assumptions would not affect the gross features observed, we hope to address these issues in 
future. 
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